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Abstract. We present details on a physical realization, in a many-body Hamiltonian system, of 
the abstract probabilistic structure recently exhibited by Gell-Mann, Sato and one of us (C.T.), 
that the nonadditive entropy Sq — k[l — Trp^]/[q~ 1] (p = density matrix; = —kTrplnp) can 
conform, for an anomalous value of q (i.e., q^ 1), to the classical thermodynamical requirement 
for the entropy to be extensive. Moreover, we find that the entropic index q provides a tool to 
characterize both universal and nonuniversal aspects in quantum phase transitions {e.g., for a 
L-sized block of the Ising ferromagnetic chain at its T = critical transverse field, we obtain 
Iknz.-^x.S' ^_^{L) / L = 3.56 ±0.03). The present results suggest a new and powerful approach to 
measure entanglement in quantum many-body systems. At the light of these results, and similar ones 
for a li = 2 Bosonic system discussed by us elsewhere, we conjecture that, for blocks of linear size L 
of a large class of Fermionic and Bosonic cZ-dimensional many-body Hamiltonians with short-range 
interaction ai T — 0, we have that the additive entropy Si{L) ^ [L''^' — l\/{d — 1) {i.e., InL for 
d — I, and V^^^ for d > \), hence it is not extensive, whereas, for anomalous values of the index 
q, we have that the nonadditive entropy Sq{L) °^ (Vt/), i.e., it is extensive. The present discussion 
neatly illustrates that entropic additivity and entropic extensivity are quite different properties, even 
if they essentially coincide in the presence of short-range correlations. 
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INTRODUCTION 

The appearance of long-range correlations in the ground state of a quantum many- 
body system, undergoing a quantum phase transition at zero temperature, is due to 
the entanglement [1]. Quantum spin chains, composed by a set of localized spins cou- 
pled through short-range exchange interaction in an external transverse magnetic field, 
capture the essence of these intriguing phenomena and have been extensively studied 
[3, 2, 4, 5, 6]. The degree of entanglement between a block of L contiguous spins 
and the rest of the chain in its ground state, as measured by the von Neumann block 
entropy Si{L) = —kTr pilnpi {pi = Trp^^ipj^ is the reduced density matrix of a L- 
sized block within a N ^ oo chain with density matrix p^v), typically saturates (i.e., 
\imi^ooSi{L) < oo) or is logarithmically unbounded (i.e., Si{L) <x InL) for large block 
size, off or at the critical point, respectively. Here we show that the nonadditive entropy 

1— TrD* 

[7, 8] S^(L) = k^-^ of the block of L spins of the ground state of quantum spin chains 
in the neighborhood of a quantum phase transition is extensive (z. e. , for L 1 , 5^ (L) L) 
for special values of ^ < 1. The additive von Neumann entropy Si{L) = —kTrp^lnpi 



is (like the additive Renyi entropy) nonextensive; indeed, liniL^ooSi (L)/L = in all 
considered cases. We present here details of the first physical realization (this as well 
as another, Bosonic, physical realization have been discussed in [9]), in a many-body 
Hamiltonian system, of the abstract mathematical examples shown in Ref. [10], that, 
for anomalous values of q, the nonadditive entropy Sq, can be extensive, as expected 
from the Clausius thermodynamical requirement for the entropy. We find that the index 
q provides a new and efficient tool to characterize different universality classes in quan- 
tum phase transitions, and to quantify entanglement [11, 12, 13] in quantum many-body 
systems, by using a nonadditive measure [14, 15, 16, 17, 18, 19, 20, 21]. 

NONEXTENSIVE STATISTICAL IMECHANICS 

The aim of statistical mechanics is to establish a direct link between the mechanical 
microscopic laws and classical thermodynamics. The most famous classical theory in 

this field has been developed by Boltzmann and Gibbs (EG) and it is considered one of 
the cornerstones of contemporary physics. The connection between micro- and macro- 
world is described by the so called BG entropy: 

w 

SBG=-kY,Pi^^Pi (1) 

i=l 

where is a positive constant, W is the number of microscopic states and {pi}i=i,...,w 
is a normalized probability distribution. One of the crucial properties of the entropy 
in the context of classical thermodynamics is extensivity, namely proportionality with 
the number of elements of the system. The BG entropy satisfies this prescription i/the 
subsystems are statistically (quasi-) independent, or typically if the correlations within 
the system are generically local. In such cases the system is called extensive. 

In general, however, the situation is not of this type and correlations may be far from 
negligible at all scales. In such cases the BG entropy may be nonextensive. Nonetheless, 
for an important class of such systems, an entropy exists which is extensive in terms of 
the microscopic probabilities [10]. The additive BG entropy can be generalized into the 
nonadditive ^-entropy [7] 

S =k ^i^, qe^ {Si = Sbg). (2) 

q-l 

This is the basis of the so called nonextensive statistical mechanics [8], which general- 
izes the BG theory. 

Additivity (for two probabilistically independent subsystems A and B) is gener- 
alized by the following pseudo-additivity: Sq{A,B)/k = Sq{A) /k + Sq{B) /k + (1 — 
q)Sq{A)Sq{B) /k^; the cases q < \ and q> \ correspond to super-additivity and suh- 
additivity, respectively. For subsystems that have special probability correlations, 
extensivity is not valid for Sbg, but may occur for Sq with a particular value of the index 
q^l. Such systems are sometimes referred to as nonextensive [10, 8]. 

A physical system may exhibit genuine quantum aspects. In particular, quantum 
correlations, quantified by the entanglement, can be present. The classical probability 



concepts are replaced by the density matrix operator p, in terms of a more general 
probability amplitude context. Therefore the quantum counterpart of the BG entropy 
in Eq. (1), which is called von Neumann entropy, is given by Si (p) = —kTr p Inp, while 
the classical ^-entropy, Eq. (2), is replaced by: 

1 -Tr/3^ 

S,{p) = k^-^. (3) 

The pseudo-additivity property is now given by 

S.iPi^Pl) _S,{pi) ^ 5,(P2) ^ ^^ Sqipi) Sq{p2) . 

k Jc k k k 

from now onk= 1 . 



XY MODEL 

In this paper we analyze a quantum system in which strong non-classical correlations 
occur between its components. We focus our investigations on a one-dimensional spin- 
1/2 ferromagnetic chain with an exchange (local) coupling and in the presence of an 
external transverse magnetic field, i.e., the quantum XY model. The Hamiltonian of the 
XY model with open boundary conditions is: 



=^ = -1 (i+r)^/^;+i + (i-r)<^JajVi+2Ad; 



(4) 



where a" (a = j, z) are the Pauli matrices of the j-th spin, N is the number of spins of 
the chain, /and A characterize, respectively, the strength of the anisotropy parameter and 
of a tras verse magnetic field along the z direction. This model for < |7| < 1 belongs 
to the Ising universality class and it actually reduces to the quantum Ising chain for 
\y\ = 1. This system undergoes a quantum phase transition at the critical point \Xc\ = 1 
in the thermodynamic limit — > oo. For 7 = it is the isotropic XX model, which is 
critical for |A| < 1 [1]. Let us stress that in the following we will solve analytically the 
ground state in the limiting case of an infinite chain, i.e. N ^00, Therefore, the coupling 
being ferromagnetic, the results will not depend on our particular choice of the boundary 
conditions. 

The entanglement in the neighborhood of the quantum phase transition has been 
recently widely investigated, thus establishing a direct connection between quantum 
information theory and condensed matter physics [2, 3, 4, 5, 6]. In particular it has been 
shown that one-site and two-site entanglement between nearest or next-to-nearest spins 
display a peak near or at the critical point [2, 3]. On the other side, the entanglement 
between a block of L contiguous spins and the rest of the chain in the ground state, 
quantified by the von Neumann entropy, presents a logarithmic divergence with L at 
criticality, while it saturates in a non-critical regime [4, 5, 6]. 

The inadequacy of the additive von Neumann entropy as a measure of the information 
content in a quantum state has been pointed out in Ref. [15]. A theoretical observation 



that the measure of quantum entanglement may not be additive has been discussed in 
Refs. [16, 14, 15, 17, 18, 19, 20, 21]. Recently, Ref. [22] suggested to abandon the a 
priori probability postulate going beyond the usual BG situation. 

Here we propose to extend the definition of the von Neumann entropy to a wider 
class of entropy measures which naturally include it, thus generalizing the notion of the 
block entanglement entropy. The block ^-entropy of a block of size L is simply defined 
as the ^-entropy, Eq. (3), of the reduced density matrix pi of the block, when the total 
chain is in the ground state. In the following we show that, contrary to the von Neumann 
entropy, there exists a q value for which Sq{pi) is extensive. This value does depend on 
the critical properties of the chain and it is consistent with the universality hypothesis. 

The XY model in Eq. (4) can be diagonalized exactly with a Jordan-Wigner transfor- 
mation, followed by a Bogoliubov rotation [23, 24, 25, 26]; this allows one to analyt- 
ically consider the thermodynamic limit N ^ The normal modes of the system are 
linear combinations of the following non-local Majorana fermions: 
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(5) 
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These operators are Hermitian and obey the anti-commutation rules {cm,c„} = 25„ 
The ground state \^g) is completely characterized by the scalar product {cmCn) 

Smn + i^iS, where 
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with real coefficients g/ given, for an infinite chain, hy 8i = j^ Jq^ d^e '^^ 



-A— lysincj 



|cos0— A— ('/sin^l ■ 

The spectrum of pi in an infinite chain in its ground state can then be exactly evalu- 
ated [6]. Indeed, the matrix pi can be written as a tensor product in terms of L uncorre- 
lated non-local Fermionic modes, which are linear combinations of the operators c„ in 
Eq. (5): pi — fi . . . fi, where f/ denotes the mixed state of mode /. The eigenvalues 
of f; are (1 ± Vz)/2, where V; is the imaginary part of the eigenvalues of the matrix r^^\ 
The entropy in Eq. (3) is then easily computed by using the pseudo-additivity, Eq. (4), 
and by noticing that the trace of ff is simply Tr = [(1 + Vz)/2]^ + [(1 - V/)/2]^. 
Notice that the required computational time scales polynomially with the block size L, 
thus allowing one to reliably analyze blocks with up to a few hundreds of spins. 



RESULTS 



We first analyze the anisotropic quantum XY model, Eq. (4) with 7^0, that has a 
critical point in Ac = 1. The block ^-entropy as a function of the block size can show 
completely different asymptotic behaviors, by varying the entropic index q. In particular, 
here we are interested in a thermodynamically relevant quantity, namely the slope. 



noted Sq, of Sq versus L. It is generically not possible to have a finite value of ^i: 
the entanglement entropy, evaluated by the von Neumann entropy, either saturates or 
diverges logarithmically in the thermodynamic limit, for respectively non-critical or 
critical spin chains [4, 5, 6]. The situation dramatically changes by using the entropy in 
Eq. (3): qualitatively it happens that, regardless the presence or absence of criticality, a 
A -dependent value of q, noted qem, exists such that, in the range I <^L<^t, {t, being the 
correlation length), is finite, whereas it vanishes (diverges) for q > qent (l < lent)- 
We note that here the nonextensivity {i.e., q ^ I) features are not due to the presence 
of say long-range interactions [27] but they are triggered only by the fully quantum 
nonlocal correlations. In Fig. 1 we show, for the critical Ising model (A = 1, 7= 1), the 
behavior of the block ^-entropy with respect to the block size: Sq{pL) becomes extensive 
(i.e., < \imi^^Sq{pi)/L<oo) for qent — 0.0828 ± 10^^ (with a corresponding entropic 
density Sq^^^ ^3.56± 0.03), thus satisfying the classical thermodynamical prescription. 
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FIGURE 1. Block ^-entropy Sq{pL) as a function of the block size L in a critical Ising chain (7 = 

1 , A = 1), for typical values of q. Only for q — qe,„ ~ 0.0828, is finite (i.e., Sq is extensive); for q < qe„t 
(q > qent) it diverges (vanishes). 

A very similar behavior is shown for non-critical Ising model, as well as for critical 
and non-critical XY models with < 7 < 1. The value of qent, for which Sq{pL) is 
extensive, is obtained maximizing numerically the linear correlation coefficient r of 
Sqipt), in the range 1 <^ L <^ ^, with respect to q, as shown in the bottom inset in Fig. 

2. Let us stress that, at precisely the critical point, ^ diverges, hence L is unrestricted 
and can run up to infinity. The index qent depends on the distance from criticality and 
it increases as A approaches (Fig. 2). It is worth stressing that our numerical results 
satisfy the duality symmetry A — > 1/A, investigated in Ref. [29]. 

We have also checked other values of 7 for the XY model and the results are very 
similar to those presented here. This fact is consistent with the universality hypothesis. 
On one hand, XY and Ising model (Ising universality class) have the same behavior as 
regards the extensivity of Sq{pi) ; in Fig. 3 we report the variation of with respect to 
A . On the other hand, for the isotropic XX model (7 = 0) in the critical region | A | < 1 we 
find qent ~ 0. 15 ± 0.01 (~ 2q^^f with ~ 0.08) for which Sqipt) becomes extensive. 
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FIGURE 2. The A-dependence of the index qent in the Ising (7=1, circle) and XY (7 = 0.75, square) 
chains. At bottom: Determination of q^m through numerical maximization of the linear correlation coef- 
ficient r of Sq{pL)- The error bars for the Ising chain are obtained considering the variation of qem when 
using the range 100 <L< 400 in the search of 5^(Pl) linear behavior Actually, at the present numerical 
level, we cannot exclude finite-size effects off criticality. 
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FIGURE 3. The A-dependence of the ^-entropic density 5^^,,, in the Ising (7= 1, circle) and XY 
(7= 0.75, square) models. For A = 1, the slopes are 3.56 and 2.63, for 7== 1 and 7= 0.75, respectively. 



Ref. [28] enables us to analytically confirm, at the critical point, our numerical results. 
The continuum limit of a (l-i-l)-dimensional critical system is a conformal field theory 
with central charge c. In this quite different context, the authors re-derive the result 
Si{Pl) ~ (c/3)lnL for a finite block of length L in an infinite critical system. To 
obtain the von Neumann entropy, they find an analytical expression for Trp|, namely 
Trpl ~ L-^/H^-i/q)_ Here, we use this expression quite differently. We impose the 



extensivity of Sq{pi) finding the value of q for which —c/6{qent — ^/qem) = 1, i.e.. 



qent = • (6) 

C 

Consequently, limi^ooS ^^—2_^{L) / L < 0°. When c increases from to infinity (see 

c 

Fig. 4), qent increases from to unity (von Neumann entropy); for c — 4 (dimension 
of physical space-time), q — 1/2; c = 26 corresponds to a 26-dimensional Bosonic 
string theory, see [30]. It is well known that for critical quantum Ising and XY models 
the central charge is equal to c = 1/2 (indeed they are in the same universality class 
and can be mapped to a free Fermionic field theory). For these models, at A = 1, the 
value of q for which Sq{pL) is extensive is given by qem — \/37 — 6 ~ 0.0828, in perfect 
agreement with our numerical results in Fig. 2. The critical isotropic XX model (7=0 
and \X\ < 1) is, instead, in another universality class, the central charge is c = 1 (free 
Bosonic field theory) and Sq{pi) is extensive for qent = \/lO — 3 ~ 0.16, as found also 
numerically. We finally notice that, in the c — >^ 00 limit, qem 1- We do not clearly 
understand the physical interpretation of this fact. However, since c in some sense plays 
the role of a dimension (see [30]), this limit could correspond to some sort of mean field 
approximation. If so, it is along a line such as this one that a mathematical justification 
could emerge for the widely spread use of BG concepts in the discussion of mean-field 
theories of spin-glasses (within the replica-trick and related approaches). Indeed, BG 
statistical mechanics is essentially based on the ergodic hypothesis. It is firmly known 
that glassy systems (e.g., spin-glasses) precisely violate ergodicity, thus leading to an 
intriguing and fundamental question. Consequently, a mathematical justification for the 
use of BG entropy and energy distribution for such complex mean-field systems would 
be more than we''^'^"^'^ 




FIGURE 4. qent versus c with the ^-entropy, >S'^(pt), being extensive, i.e., liniL-^^S {Pl)/L < «>. 

c 

When c increases from to infinity, qem increases from to unity (von Neumann entropy); for c = 4, 
q = 1/2 and for c 3> 1, see Ref. [30]. Inset: for the critical quantum Ising and XY models c = 1/2 and 
qent = — 6 ~ 0.0828, while for the critical isotropic XX model c = 1 and qem = a/TO — 3 ~ 0. 16. 



It is worth to mention that the Renyi entropy of a block of critical XX spin chains has 
been derived analytically in Ref. [31]. Since the Renyi entropy is simply connected to 
the entropy Sq, it is possible to re-derive qem for the critical XX model also from that 
analytical expression. 

CONCLUDING REMARKS 

Summarizing, we have presented: (i) Details on the first physical realization (in a 1/2- 
spin d = I quantum system), in a many-body Hamiltonian system, of the abstract 
probabilistic structure shown in Ref. [10], that Sq conforms, for a special value of 
q, to the classical thermodynamical requirement for the entropy to be extensive (the 
second physical realization, in a = 2 Bosonic system, can be seen in [9]); (ii) A 
new connection, Eq. (6), between nonextensive statistical mechanical concepts and BG 
statistical mechanics at criticality (see [32] for another such analytical connection); (iii) 
A novel and simple manner to characterize entanglement through the pair (qenti Sqem)- 

Let us point out also that the reduction of the pure ground state of the full chain (at 
r = 0) to a finite block of L spins results in a mixed state with quantum fluctuations. 
A mapping of this subsystem within a zero temperature XX infinite chain to a finite 
system which is thermalized at some finite temperature has been recently exhibited 
[33], thus defining an L-dependent effective temperature of the block. The use of a 
non-Boltzmannian distribution {e.g., the one emerging within nonextensive statistical 
mechanics) might enable defining an effective temperature which would not depend on 
L, as physically desirable. Indeed, this approach has been successfully implemented for 
e — e^ collision experiments [34]. 

Finally, let us emphasize the difference between additivity and extensivity for the en- 
tropy. Additivity only depends on the mathematical definition of the entropy; therefore, 
S\ is additive, whereas Sq (<7 7^ 1) is nonadditive. Extensivity is more subtle, since it 
also depends on the specific system. The T = block entropies of the present 1 /2-spin 
d —\ quantum system at criticality are given by Si{L) oc InL {i.e., nonextensive), and 
"^[Vq+c^-sJ/c^^) 0= L {i.e., extensive). It is known (see [35] and references therein) that, 
for (i-dimensional Bosonic systems (e.g., a black hole [36]), 5i follows the area law, 
i.e., S\{L) oc L^^^ {i.e., nonextensive). A logarithmic behavior for J = 1, and the area 
law for d > 1 can be unified through Si(L) oc [L^-i - l]/(j - 1) = Inz-^L [37] {i.e., 
nonextensive), which would correspond to a large class (yet not completely identified) 
of fully entangled quantum systems. For all these systems, one could expect that a value 
of q exists such that Sq{L) oc Lf^ {i.e., extensive). In addition to the present example, a 
d = 2 Bosonic system has been shown [9] to satisfy this conjecture. 

ACKNOWLEDGIVIENTS 

The present work has benefited from fruitful comments by R. Fazio, V. Giovannetti, D. 

Patane, A. Pluchino, A. Rapisarda and D. Rossini. One of us (C.T.) also acknowledges 
enlightening conversation with G. 't Hooft on the black-hole area law. This work was 
partially supported by Centro di Ricerca Matematica E. De Giorgi of the Scuola Normale 



Superiore, and by Pronex/MCT, CNPq and Faperj (Brazilian Agencies). In the 100th 
anniversary of his birthday, we dedicate this work to the memory of Ettore Majorana. 

REFERENCES 

1. S. Sachdev, Quantum Phase Transitions (Cambridge Univ. Press, Cambridge, 2000). 

2. A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 416, 608 (2002). 

3. T.J. Osborne and M.A. Nielsen, Phys. Rev. A 66, 0321 10 (2002). 

4. G. Vidal, J.I. Latorre, E. Rico and A. Kitaev, Phys. Rev. Lett. 90, 227902 (2003). 

5. A.R. Its, B.-Q. Jin and V.E. Korepin, /. Phys. A 38, 2975 (2005). 

6. J.I. Latorre, E. Rico and G. Vidal, Quant. Inf. and Comp. 4, 48 (2004). 

7. C. TsalUs, J. Stat. Phys. 52, 479 (1988). 

8. J.R Boon and C. TsalUs, eds., Nonextensive Statistical Mechanics: New Trends, New Perspectives, 
Europhysics News 36 (6) (2005). 

9. F. Caruso and C. Tsallis, preprint (2007). 

10. C. Tsallis, M. Gell-Mann and Y. Sato, Proc. Nat. Acad. Sci. 102, 15377 (2005). 

1 1 . M. Nielsen and I. Chuang, Quantum Computation and Quantum Information (Cambridge Univ. Press, 
Cambridge, 2000). 

12. V. Vedral and M.B. Plenio, Phys. Rev A 57, 1619 (1998). 

13. P Horodecki, R. Horodecki and M. Horodecki, Acta Phys. Slov 48, 141 (1998). 

14. K. Zyczkowski and I. Bengtsson, Ann. Phys. 295, 1 15 (2002). 

15. C. Brukner and A. Zeilinger, Phys. Rev A 63, 0221 13 (2001). 

16. F. Giraldi and R Grigolini, Phys. Rev A 64, 032310 (2001). 

17. C. Tsallis, S. Lloyd and M. Baranger, Phys. Rev A 63, 042104 (2001). 

18. S. Abe, Phys. Rev A 65, 052323 (2002). 

19. S. Abe and A.K. Rajagopal, Phys. Rev A 60, 3461 (1999). 

20. A.K. Rajagopal and R.W. Rendell, Phys. Rev A 72, 022322 (2005). 

21. M.B. Plenio and S. Virmani, Quant. Inf Comp. 7, 1 (2007). 

22. A. Popescu, A.J. Short and A. Winter, Nature Physics 2, 754 (2006). 

23. E. Lieb, T. Schultz and D. Mattis, Ann. Phys. 16, 407 (1961). 

24. R Pfeuty, Ann. Phys. 57, 79 (1970). 

25. E. Barouch, B.M. McCoy and M. Dresden, Phys. Rev A 2, 1075 (1970). 

26. E. Barouch and B.M. McCoy, Phys. Rev A 3, 786 (1971). 

27. C. Anteneodo and C. Tsallis, Phys. Rev Lett. 80, 5313 (1998). 

28. R Calabrese and J. Cardy, J. Stat. Mech.: Theor Exp. P06002 (2004). 

29. R. Savit, Rev Mod. Phys. 52, 453 (1980). 

30. P. Ginsparg and G. Moore, Lectures on 2D String Theory (Cambridge Univ. Press, Cambridge, 1993). 

31. B.-Q. Jin and V. E. Korepin, J. Stat. Phys. 116, 79 (2004). 

32. A. Robledo, Molecular Physics 103, 3025 (2005). 

33. V. Eisler, O. Legeza and Z. Racz, / Stat. Mech.: Theor. Exp. PI 1013 (2006). 

34. I. Bediaga, E.M.F Curado and J. M. de Miranda, Physica A 2'&6, 156 (2000). 

35. T. Barthel, M.-C. Chung and U. Schollwock, Phys. Rev A 74, 022329 (2006). 

36. M. Srednicki, Phys. Rev Lett. 71, 666 (1993). 

37. The function In^ x = ' ^^^^ {x>0; q& ^) naturally emerges within nonextensive statistical mechan- 
ics. 



